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Abstract 

In  this  study,  a  method  for  determining  the  initial  flaw  size  in  a  particulate  composite  material  is 
developed.  The  results  of  analyses  indicate  that  the  initial  flaw  length  and  the  critical  flaw  size  follow  the 
second  asymptotic  distribution  of  maximum  value.  It  also  indicates  that  on  the  first  approximation,  it  can 
be  assumed  that  the  initial  flaw  size  and  the  critical  flaw  size  is  independent  of  specimen  thickness. 

Introduction 

An  important  engineering  problem  in  structural  design  is  evaluating  structural  integrity 
and  reliability.  It  is  well  known  that  structural  strength  may  be  degraded  during  its  design 
life  due  to  mechanical  or  chemical  aging,  or  a  combination  of  these  two  aging 
mechanisms.  Depending  on  the  structural  design,  material  type,  service  loading,  and 
environmental  condition,  the  cause  and  degree  of  strength  degradation  due  to  the  different 
aging  mechanisms  differs.  One  of  the  common  causes  of  strength  degradation  is  the  result 
of  crack  development  in  the  structure.  When  cracks  occur,  the  effects  of  crack  sizes  and 
the  rate  of  growth  on  the  fracture  resistance  of  the  material  need  to  be  investigated. 

Crack  propagation  in  particulate  composite  materials  has  been  an  important  subject  of 
research  for  many  years.  It  has  been  shown  in  recent  literature  that  the  crack  propagation 
rate  in  particulate  composites  is  a  power  function  of  the  Mode  I  stress  intensity  factor  ( 
Liu,  1990;  Liu  1990;  Liu  and  Yang,  1994).  In  order  to  predict  the  crack  growth  rate,  it  is 
necessary  to  determine  the  initial  crack  size. 

It  is  well  known  in  the  aerospace  industry  that  the  initial  flaw  size  in  metals  and 
superalloys  is  too  small  to  be  detected  by  any  NDI  techniques.  Consequently,  the  initial 
flaw  size  in  metals  has  been  determined  using  experimental  results,  such  as  fractographic 
data  or  S-N  data  (Yang  et.  al,1995  and  Yang  et.  al,1997).  From  the  experimental  S-N 
data,  one  can  determine  the  terminal  crack  size  (critical  crack  size)  at  the  time  of  failure. 
The  initial  flaw  size  is  then  computed  from  the  terminal  crack  size  by  conducting  the 
crack  growth  analysis  backwards. 
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To  date,  the  concept  of  predicting  the  crack  growth  rate  based  on  the  initial  flaw  size  and 
the  determination  of  the  initial  flaw  size  in  highly  filled  polymeric  materials  has  not  been 
studied  yet.  Therefore,  it  is  the  objective  of  this  study  to  investigate  the  concept  of 
predicting  the  crack  growth  rate  based  on  the  initial  flaw  size  and  to  determine  the  initial 
flaw  size  existing  in  a  particulate  composite  material. 

In  this  study,  the  equivalent  initial  flaw  size  (EIFS)  in  a  particulate  composite  material, 
containing  hard  particles  embedded  in  a  rubber  matrix,  was  determined  using  constant 
strain  rate  crack  propagation  test  data.  Uniaxial  tensile  specimens  with  and  without  pre¬ 
crack  were  tested  at  a  constant  strain  rate  of  0.067  in/in/min.  Four  different  specimen 
thicknesses  (0.2  in.,  0.5  in.,  1.0  in,,  and  1.5  in.)  were  considered.  The  effect  of  specimen 
thickness  on  the  EIFS  was  investigated  and  the  results  are  discussed. 

Analytical  Analysis 

To  determine  the  EIFS,  the  following  information  is  needed:  (1)  crack  growth  rate 
parameters,  (2)  critical  stress  intensity  KIC  and  threshold  stress  intensity  factor  K,h  under 
which  crack  will  not  grow,  and  (3)  time  to  failure  data  under  constant  strain  rate.  Crack 
growth  rate  parameters  as  well  as  KIC  and  K,h  are  determined  experimentally  using  pre¬ 
flawed  specimens.Time  to  failure  data  are  also  obtained  experimentally  using  specimens 
without  a  pre-flaw. 

For  pre-flawed  specimens,  the  stress  intensity  factor  K,is  given  by 

Kj=  a  (n  a)1/2  f  (a/w)  (1) 

In  which  ,cr  is  the  applied  stress,  f(a/w)  is  the  geometric  correction  factor,  a  is  the  crack 
length,  and  w  is  the  width  of  the  specimen.  The  functional  relationship  between  f(a/w) 
and  a/w  is  shown  below. 

f(a/w)  =  0.5854(a/w)3+1.099(a/w)2+0.8672(a/w)+1.049  (2) 

For  a  specimen  subject  to  a  constant  strain  rate,  the  stress  intensity  factor  K,  reaches  the 
critical  stress  intensity  factor  KIC  at  the  instant  of  fracture,  and  the  corresponding  flaw 
size  is  denoted  by  ac,  referred  to  as  the  critical  flaw  size  or  the  terminal  flaw  size.  It 
follows  from  Eq.  (1)  that 

K[C=  a  c  (7rac)1/2  f  (ac/w)  (3) 

Where  a  c  is  the  critical  stress  at  fracture. 

The  crack  growth  rate  da/dt  has  been  shown  to  be  a  power  function  of  the  stress  intensity 
factor  K„  i.e., 


da/dt  =  Q  K,m 


(4) 


in  which  m  and  Q  are  crack  growth  rate  parameters. 


When  a  specimen  without  pre-flaw  is  subjected  to  a  constant  rate,  the  entire  loading 
history  and  hence  the  stress  history  a  =  cr(t)  can  be  measured,  including  the  critical  stress 
a  c  at  the  time  of  fracture,  tc.  For  a  given  critical  stress  intensity  factor  KIC  (material 
constant),  the  critical  flaw  size  ac  can  be  computed  from  Eq.  (3). Consequently,  the  initial 
flaw  size  a0  at  t=0  can  be  obtained  by  integrating  Eq.  (4),  based  on  the  terminal  condition 
(ac,  tc)  and  the  stress  history  <r(t). 

Experimental  Analysis 

Constant  strain  rate  tests  were  conducted  on  specimens  with  and  without  pre-flaws  at  a 
strain  rate  of  0.067  in/in/min.The  critical  stress  ctc  and  the  time  to  failure  tc  were 
determined  from  the  specimen  without  pre-flaw.  The  crack  growth  parameters  m  and  Q 
were  determined  form  the  specimens  with  pre-flaw.  The  results  are:  m  =  2.084  and  Q  = 
9.3325x1  O'7  in  which  the  units  are  force  in  pound,  length  in  inch,  and  time  in  minute. 
Further,  the  critical  stress  intensity  factor  and  the  threshold  stress  intensity  factor  are  78.3 
psi  (in)I/2  and  52  psi  (in)I/2,  respectively. 

Statistical  Distribution  of  Equivalent  Initial  Flaw  Size  and  Critical  Flaw  Size 

The  results  of  the  analysis  show  that  the  equivalent  initial  flaw  size  a0  as  well  as  the 
critical  flaw  size  ac  vary  from  specimen  to  specimen.  Hence,  the  statistical  distribution  of 
these  quantities  should  be  determined.  In  this  study,  four  statistical  distribution  functions 
,  (1)  normal  distribution,  (2)  two  parameter  Lognormal  distribution,  (3)  two  parameter 
Weibull  distribution  and  (4)  second  asymptotic  distribution  of  maximum  values,  were 
considered.  The  goodness  of  fit  for  different  distributions  has  been  conducted  using  the 
Kolomogorov-Smirov  test.  The  results  indicate  that  the  second  asymptotic  distribution  of 
the  maximum  value  has  the  best  fit  for  the  distribution  of  a0  and  ac. 

Results 

Based  on  the  constant  strain  rate  test  results  for  specimens  with  and  without  preflaws,  a 
method  for  determining  the  equivalent  initial  flaw  size  existing  in  the  particulate 
composite  material  has  been  presented.  The  method  is  based  on  a  backward  crack  growth 
analysis  starting  from  the  given  terminal  flaw  size  to  determine  the  equivalent  initial  flaw 
size. 

Equivalent  initial  flaw  sizes  were  computed  for  the  cases  in  which  the  effect  of  Kth  is 
ignored.Typical  plots  of  statistical  distributions  of  a0  are  shown  in  Figs.l-4.For  a 
comparison  purpose,  experimental  data,  shown  as  circles,  are  also  included  in  these 
figures.  It  is  seen  that  the  second  asymptotic  distribution  of  maximum  value  fits  the 
experimental  data  the  best  that  is  consistent  with  the  results  of  the  goodness  of  fit 
analyses.The  results  of  statistical  analyses  of  ac  are  similar  to  that  of  a0. 


Second  Asymptotic  Distribution  of  Maximum  Value 
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The  effect  of  specimen  thickness  on  the  initial  flaw  size  and  the  critical  flaw  size  were 
investigated.The  results  of  analyses  show  that  for  specimen  thicknesses  equal  to  0.2  in., 
0.5  in.,  and  1.0  in.,  the  variations  among  a0  and  ac  are  small  and  the  average  values  a0  and 
ac  are  0.126  in.  and  0.132  in.  For  specimen  thickness  equal  to  1.5  in.  the  average  values  of 
a0  and  ac  are  0.144  in.  and  0.158  in.,  and  they  are  14%  and  20%  higher  than  that  when  the 
specimen  thickness  is  equal  to  or  less  than  1.0  in.  The  increase  in  a0  and  ac  in  thicker 
specimen  is  probably  due  to  the  size  effect  of  this  material.  Unlike  the  brittle  metallic 
materials,  this  high  toughness  particulate  composite  material  shows  higher  strength  when 
the  size  of  the  specimen  is  increased.  The  size  effect  of  this  material  is  a  subject  of  further 
research.  Considering  the  highly  nonhomogeneous  nature  of  the  highly  filled  particulate 
composite  and  for  the  engineering  application  purpose,  it  is  reasonable  to  assume  that  a0 
and  ac  are  independent  of  the  specimen  thickness.  Under  this  assumption,  the  average 
value  and  the  coefficient  of  variation  of  a0  (ac)  are  0.131  in.  (0.146  in.)  and  0.0703  (0.054 
),  respectively. 

Conclusion 

In  this  study,  a  method  is  developed  to  predict  the  initial  crack  length  and  its  statistical 
distribution  function  of  a  particulate  composite  material  subjected  to  a  constant  strain  rate 
loading  condition.  The  results  of  analyses  indicate  that  the  equivalent  initial  flaw  size  and 
the  critical  flaw  size  follow  the  second  asymptotic  distribution  of  maximum  value  .It  also 
indicates  that,  on  the  first  approximation  and  for  the  engineering  application  purpose,  it 
can  be  assumed  that  the  initial  crack  length  is  independent  of  specimen  thickness. 
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